Linearons: highly non-instantaneous solitons in liquid-core photonic crystal fibers 
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The nonlinear propagation of light pulses in liquid-filled photonic crystal fibers is considered. Due 
to the slow reorientational nonlinearity of some molecular liquids, the nonlinear modes propagating 
inside such structures can be approximated, for pulse durations much shorter than the molecular 
relaxation time, by temporally highly-nonlocal solitons, analytical solutions of a linear Schrodinger 
equation. The physical relevance of these novel solitary structures, which may have a broad range 
of applications, is discussed and supported by detailed numerical simulations. 



Introduction — Highly nonlocal solitons have been 
originally introduced as an "accessible" toy-model for de- 
scribing self-trapped optical beams [l[; but this was fol- 
lowed by experimental demonstrations , also driven by 
early works on plasma physics, Bose-Einstein condensa- 
tion and dissipative systems |3l4lCj{. which unveiled the 
fundamental role of nonlocality in spatially self-trapped 
waves. Indeed, nonlocality allows stabilization with re- 
spect to collapse and the existence of a rich class of 
propagation-invariant waves lllll2j : in addition, applica- 
tions such as light-steering and all-optical logic gates [l3| 
and in soft-matter and thermal liquids jl4 - |l6l ] have been 
demonstrated. Recently, insights from spatially nonlocal 
nonlinear waves also emerged in the temporal domain 
17 1. However, the relevance of temporal nonlocality is 
largely limited by the unavoidable instantaneous Kerr 
effect, as, e.g., for silica glass in fiber optics, where "non- 
local" Raman-like terms, although leading to important 
consequences as the Raman self-frequency shift (RSFS) 
of solitons [li| , can be considered as a small perturbation. 
Recent fabrication advances, however, open up innova- 
tive perspectives. Indeed, micro-structured photonic- 
crystal fibers (PCFs) may be fabricated with a central 
hole filled by a material displaying non-instantaneous re- 
sponse as, e. g., molecular liquids with large reorienta- 
tional effects [19j, with instantaneous nonlinearities act- 
ing as small perturbations: exactly the opposite situation 
of the silica glass. 

In this Letter we show that a completely novel class 
of solitary wave exists in these new fibers, which are de- 
scribed by an essentially linear model. These linearons 
are shown to display truly remarkable properties, such as 
robustness with respect to noise perturbations and colli- 
sions, and emission of phase-matched dispersive waves - 
phenomena that are typically associated to purely non- 
linear temporal waves only. In addition they are shown 
to be not affected by the RSFS if the response time of 
the molecular reorientation is slow enough. As a result, 
these highly noninstantaneous fibers can be used in a new 
variety of optical devices and applications, also including 
quantum information processing. 

Liquids with reorientational nonlinearities — Suit- 
able liquids with large reorientational nonlinearities and 



comparably small Kerr contributions are those com- 
posed of small molecules. Good candidates are sol- 
vents typically used for spectroscopy. These liquids, 
which are available in extremely high purity, show 
large transparency windows, ranging from 500 nm to 
3 fim [2^], high optical damage threshold and compa- 
rably long reorientational response times. Experimen- 
tal implementation requires the liquid to possess cer- 
tain physical properties in addition to the nonlinear op- 
tical ones. An in-depth survey shows that best po- 
tential candidates are nitrobenzene (C6H5NO2), toluene 
(CyHg) and carbon disulfide (CS2), which all pos- 
sess small cigar-shaped molecules with comparable low 
Kerr nonlinearity (C 6 H 5 N0 2 : 671xlO- 16 cm 2 /W, C 7 H 8 : 
168xlO- 16 cm 2 /W, CS 2 : 340x 10" 16 cm 2 /W; all at 1064 
nm [2l|) and refractive indices larger than silica. From 
a practical point of view, both nitrobenzene and toluene 
have relatively low vapor pressure at room temperature 
and hence are easy to handle. CS2 , in contrast, is volatile, 
but has the best transparency properties due to its sim- 
pler molecular composition. To illustrate our findings 
we have performed calculations for CS2, which is a com- 
mon reference material in spectroscopy and thus all of its 
optical and physical properties are extremely well docu- 
mented (22)]. However, we have to stress that our results 
are of general character and are not restricted to one 
specific liquid. 

Linearon equation — Our starting point is the GNLSE 
written in dimcnsionless form: 



id z A + -sd?A + A 



+ OO 



R{t-t')\A(t')\ 2 dt' = 0, (1) 



where R(t) = (1/T) exp(— t/T)Q(t) is the response func- 
tion of the CS2 reorientational nonlinearity, T = 1.68 
psec is the medium response time (22[, 0(f) is the Heavi- 
side function that ensures causality, A(z, t) is the electric 
field envelope, s = +1 and s = — 1 denote anomalous and 
normal group velocity dispersion (GVD) respectively, z is 
the dimensionless propagation coordinate along the fiber, 
and t is the dimensionless time coordinate. The response 
function R(t) must be normalized to one: J R(t)dt = 1. 
We further define the moments of the response function: 
(T(™)) eee / t m R(t)dt. 
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In the well-known limit (valid for materials such as sil- 
ica, for instance, see [IH) of first moment much shorter 
than the pulse duration to, i.e. T = T 1 - 1 ^ -c to, one 
can use the reversibility property of the convolution in- 
tegral R(t - t')\A(t')\dt' = R(t')\A(t - t')\dt' 
to expand the envelope in a Taylor series: \A(t — t')\ 2 ~ 
\A(t)\ 2 - t'd t \A(t)\ 2 + \t"d*\A(t)\ 2 + .... By inserting 
this expression into Eq . dJ) , we obtain the well-known 
Raman-NLSE model 0: id z A + \sd^A + \A\ 2 A - 
TAd t \A\ 2 = 0, where the last term models the dynam- 
ics of long pulses subject to the Raman effect. In sil- 
ica fibers this approximation is quite good for picosec- 
ond pulses, since T ~ 2 fsec, [18(. A completely dif- 
ferent scenario occurs in the opposite limit, i.e. when 
T > t - In CS 2 liquid, T ~ 1.68 psec [§2- In 
that case, one expands in a Taylor series the response 
function (which is varying slowly with respect to A): 
R(t-t') ~ R{t)-t'd t R{t) + \t"d 2 R{t) + .... Maintaining 
only the zero-th order term, Eq. (TTJ) becomes: 



found by solving the implicit equation 



id z A + -sd 2 A 



£R(t)A = 0, 



(2) 



where £ = /_ °° |^4(i)| 2 di is the total number of photons 
launched into the fiber. In the highly non-instantaneous 
limit, the GNLSE (pj) is equivalent to a linear Schrodinger 
equation with a time-dependent potential R{t). In this 
regime, the response time is so long with respect to the 
pulse duration, that the system 'remembers' the total 
energy injected into the fiber, while the pulse is prop- 
agating. Consequently, the pulse feels the exponentially 
decaying response of the liquid as it were a purely linear 
potential, whose depth is determined by the energy £. 

Linearons states — Localized solutions of Eq. ^ can 
be found by solving the equation for t < and t > 0, 
and then by imposing continuity of the envelope and its 
derivative at t = 0. These continuity conditions provide 
the constraint equation for the wavenumber (3 > 0. The 
solution of @ with the exponential potential is known 
from standard quantum mechanics, being the exponential 
well among the solvable potentials. Soliton solutions of 
Eq. @ arc found by posing A(z,t) — a(t) exp (if3z), 
which leads to (for t > 0) 



^d 2 a+^eM-t/T)a = /3a. 



(3) 



For (t < 0) the exponentially localized solution is a = 
A/"exp(v2/3£) where Af is a normalization constant, de- 
termined by the total energy £. Among the possible so- 
lutions of ([3]) for t > 0, we choose those vanishing as 
t — > oo, which leads to 



Af 



J v (V8£Te 



-t/2T\ 



J V (V8£T) 



(4) 



with v = y/8(3T. For large t, a decays as exp(—vt/2T). 
The allowed eigenvalues of /3 for a specified value of £ are 



o^i (; 



-2St) = 0, 



(5) 



where o-fiG a; z ) = S^Lo z " /[ n K a )n] is the confluent hy- 
pergeometric limit function, and (a) n = T(a + n)/T(a) 
is the rising factorial. Eq. ([5]), which provides the dis- 
persion relation (3{£,T), can only be solved numerically. 
The normalization constant Af, which cannot be written 
explicitly, is found by requiring that the total soliton en- 
ergy to equal £. Solutions can be classified by counting 
the number of nodes (m). We shall call the 'fundamen- 
tal' state the solution with m = 0, which has the largest 
j3 and the shortest pulse duration. For small values of £ , 
Eq. can be expanded to give (3 ~ £ 2 /2, which shows 
that /3 = £ = is also a solution. For large values of £, 
we have found that, for the fundamental state of largest 
f3, one has the quasi-linear behavior j3 ~ (2/ir) 2 £ /T . An- 
alytical solutions for moving solitons can also be found 
by taking the solutions of Eq. §3§ and making the sub- 
stitution 2/3 — > 2/3 — v 2 , and by multiplying the fields 
by a factor e~ lvt , where v is an extra parameter associ- 
ated with the soliton velocity. In this case, we have an 
extra constraint j3 > v 2 /2 that must be satisfied for a 
real solution to exist. For \t\ large, the field decays as 
a ~ exp(— \/2^|t|). This means that the original approx- 
imation used to derive Eq. ([2J from Eq. (TTJ is verified a 
posteriori only if /3 3> [2T 2 ] -1 , or by using the large en- 
ergy limit £ ^> 7r 2 /(8T), which is the physical condition 
for our simplified highly-nonlocal model to be valid. 

Figure H{a) shows the profiles of some analytical solu- 
tions of Eq. (j2J|. Figures QJb,c) show the dependence of 
/3 on £ and T for fixed T and £ respectively, as calculated 
by solving Eq. ((SJ) numerically. Approximate expressions 
in terms of £ and T are given below. 

Raman self-frequency shift of linearons — The analyt- 
ical solutions of the reduced model ^ calculated above, 
being linear, are stable against small noise fluctuations 
in their profile. To prove this, in Fig. [2ja) we show 
as an example the propagation of the higher-order 
soliton shown in Fig. [lja) [blue line], when perturbing 
it with 5% noise. If the same analytical solution is 
propagated inside the full model (TTJ), we observe that 
the pulse is subject to a Raman shift, which is not 
described by the simplified model of Eq. @. We 
have developed a general theory of RSFS of solitons 
for Eq. (Qj). We rewrite Eq. (JT]) as id z A + \sdfA + 
£R(t)A + [J R(t - t')\A(t')\ 2 dt' - £R(t)] A = 0, 
and treat the term under square brackets as a 
perturbation of the stationary states found in the 
previous section. One can prove that the central 
frequency of the soliton shifts according to Q,(z) = 
-(z/£)Jdt\A(t)\ 2 d t [J R{t - t')\A(t')\ 2 dt' - £R(t)] ~ 
c As zJdt[d 2 \A(t)\ 2 ]R{t), where c AS = 
Jt'\A(t')\ 2 dt'/ J\A(t')\ 2 dt' is the asymmetry coefficient 
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FIG. 1: (Color online) (a) Plots of localized solutions of Eq. 
© for T = 10, £ = 30, for (3 = 2.33614 (m = 0), 1.82551 
(m = 1), 1.45665 (m = 2) and 1.16691 (m = 3), calculated by 
using Eq. ©. (b) Af, cas, <?, P as functions of T, for £ — 30. 
(c) Plot of P(T) for fixed £ = 10. (d) Plot of /3(£) for fixed 
T= 10. 



of the localized state. The final result is (x = V8£T): 



n 



J„(x) 2 T JO e 



■ t / T ^j2( a . e -t/2T )dt 



32 f- 1 



(6) 

which shows that the RSFS of these solitons depends 
on the asymmetry of their temporal profiles. The last 
expression in ([6]) is valid for large £ and T (where it is 
possible to show M 2 = 2£/tt 2 T and cas — T/n 2 ), and 
illustrates how, as the degree of nonlocality T grows, the 
Raman shift is inhibited, while also growing with the 
energy. The rate of RSFS g = Cl/z and the asymmetry 
coefficient cas of the stationary states are shown in Fig. [5] 
in terms of T and £ for the fundamental soliton solutions. 

Dispersive resonant radiation — In the presence of 
higher-order dispersion terms in Eq. ([2]), one can show 
that linearons can resonantly emit dispersive radiation 
at well-defined frequencies, analogously to what occurs 
for Schrodinger solitons in solid-core PCFs 23]. By 
substituting A(z,t) = [F(t) + f{z 1 t)]e l ^ z into Eq. ©, 
where F(t) is the linearon profile and / is the small 
dispersive radiation amplitude, assuming that the re- 
sponse time T is large in comparison with the linearon 
duration, and keeping only the first order term in /, 

we have: id z - (3 + D(id t ) + £R(t) f = -D R (id t )F, 

where D{id t ) = (l/2)d 2 + iadf and £>h = iadf, and 
a is the third-order dispersion coefficient, the only one 
that we include here. At phase-matching, radiation 
and soliton have the same wavenumber /3 and f(co) — 
S{lu)/[D(lu)+£R(lu)- P], where S(lu) is the Fourier 
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FIG. 2: (Color online) (a) Propagation of 1-soliton in full eq., 
T = 10, E = 30, P = 2.33614. (b) Same as (a) but for 
T = 250, £ = 30, P = 0.109846, showing reduction of RSFS- 
induced acceleration, (c) Various final spectra showing re- 
duction of RSFS (£ = 32, T = 10, 50, 100, 250). (d) Resonant 
radiation in reduced model for a — 0.05 (for m — 4 linearon). 
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FIG. 3: (Color online) (a,b) Temporal oscillations of linearons 
due to a small instantaneous Kerr nonlinearity [r — in (a) 
and r — 1/6 in (b), T = 10, input pulse A = iVsech (t) 
with N = 5], in the reduced model of Eq. ©. (c) GVD 
of the fundamental core mode (HEn-mode) of a CS2-filled 
PCF. The structure consists of the liquid CS2 core and one 
additional ring of air holes (see the inset), with a = 1.6 /im, 
b = 1.2 fim. Each curve refers to a structure of a different 
hole spacing A expressed in fj,m. 
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transform of the source term —Dn(idt)F . This yields 
the phase-matching condition D(lu) +£R(uj) ~ f3, which 
determines the resonant frequency wr. The energy- 
dependent part, although quite small for T — > oo, is 
an extra contribution to the resonant condition that is 
unique for linearons, and allows the frequency position 
of the emitted radiation to be tuned by adjusting the 
total input pulse energy injected into the fiber. 

Kerr nonlinearity — An important issue concerns the 
effect of a residual instantaneous Kerr nonlinearity on 
the linear dynamics. This may be due to the Kerr effect 
of the fiber cladding, or of the liquid itself. As shown 
in panels [3^a,b) for a sech-like input pulse (i.e. an input 
that is not matched with the profile of the fundamental 
linearon), in the absence of the Kerr effect [Fig. EJa)], 
the propagation largely resemble a standard fundamen- 
tal soliton. Conversely, increasing the contribution of the 
Kerr effect, parameterized by a coefficient r which is the 
relative importance between the instantaneous and the 
non-instantaneous part of the nonlinearity (typically of 
the order of ~ 0.1 in liquids) [Fig. [3(b)], induces a char- 
acteristic oscilation of the central position of the pulse 
in the temporal and in the spectral domain. Such an ef- 
fect is due to the coupling between higher order linearons 
that get excited in the fiber, induced by the perturbation 
introduced by the Kerr nonlinearity. The details of such 
a perturbation theory will be reported elsewhere. 

Linearon-fiber — In Fig. [3jc) we show the calculated 
GVD of a series of CS2-filled core silica PCFs (with a sin- 
gle ring of holes arranged in a triangular cladding lattice 
as shown in the inset) for different values of the pitch A. 
Other parameters are given in the caption. The PCF in- 
troduces tremendous flexibility in the engineering of the 
position of the zero-GVD point, which can be shifted over 
a large range of wavelengths. This will allow optimal ex- 
perimental excitation of linearons in the fiber, and the 
formation of linearon-induccd supercontinua. 

Conclusions — We predict the existence of a novel 
class of temporally localized waves propagating inside mi- 
crostructured fibers with a central core filled by nonlin- 
ear liquids with a slow reorientational nonlinearity. Sur- 
prisingly, these waves (which we dubb linearons) behave 
very much like solitons - albeit being the solution of a 
linearized equation - and are sustained by a strongly non- 
local temporal response, due to a pronounced Raman-like 
effect induced by reorientational nonlinearity of cigar- 
shaped molecules. By borrowing concepts from nonlo- 
cal spatial solitons, we find that these objects are very 
robust with respect to noise, caused by for example am- 
plified spontaneous emission. Linearons may thus sup- 
port the development of novel classes of lightwave fiber 
systems and novel soliton based lasers, as well as open- 
ing up a new route towards quantum solitons and mul- 



tidimensional solitary waves sustained by a highly non- 
instantaneous nonlinearity. In addition, by exploiting the 
nonlinear coupling between linearons, we believe that it 
will be possible to control and enhance supercontinuum 
generation in novel and unexpected ways. 
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